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Abstract 

We consider localized states in a discrete bistable Allen-Cahn equation. This model equation combines 
bistability and local cell-to-cell coupling in the simplest possible way. The existence of stable localized states 
is made possible by pinning to the underlying lattice; they do not exist in the equivalent continuum equation. 

In particular we address the existence of 'isolas': closed curves of solutions in the bifurcation diagram. 
Isolas appear for some non-periodic boundary conditions in one spatial dimension but seem to appear 
generically in two dimensions. We point out how features of the bifurcation diagram in ID help to explain 
some (unintuitive) features of the bifurcation diagram in 2D. 
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1. Introduction 

The general term 'pattern formation' refers to the spontaneous appearance of structure in a dissipa- 
tive dynamical system in a continuum or spatially discrete medium, when a spatially homogeneous state 
undergoes an instability of some kind [12] [I7j [22]. In many physical situations it is appropriate to use 
mathematical models for which the spatial dynamics are modelled discretely, rather than as a continuum. 
Such lattice models arise in many physical applications: in particular we mention recent work in nonlinear 
optics [28] [29] [8] and mathematical modelling of processes in developmental biology [10] [2TJ [26]. In many 
applications the existence of domain-filling (almost) regular solutions exhibiting spatial periodicity is of in- 
terest. Mathematical work on periodic patterns in regular lattices has also been the subject of considerable 
recent interest [25] [2] . 

In addition, there has been a renewed interest in the study of localized states in nonlinear dissipative 
systems in recent years, motivated by a wealth of experimental and numerical results both in continuum 
systems [HJ[4j[5] anc ^ s P a ti a lly discrete models [1 . Such localized states have been referred to by a number 
of different names, including 'light bullets' and 'dissipative solitons'. 

In continuum systems one recurring mechanism for the generation of localized states is the existence 
of a subcritical Turing (pattern- forming) instability. In the simplest case, two branches of small- amplitude 
localized states emerge at the Turing instability, also in this context referred to as a Hamiltonian-Hopf 
bifurcation [18] . Away from the initial instability they develop into a complicated structure of fully nonlinear 
localized states which exist on a pair of intertwining curves [27 ; a complete asymptotic description of the 
bifurcation structure involves consideration of exponentially-small effects [19] [7] . In an infinite domain this 
'homoclinic snaking' behaviour continues, and the width of the localized states increases without bound. 
In a finite domain the branches of localized states terminate on the branches of spatially periodic patterns 
which also originate at the initial Turing instability [3] [13] . There is now a pretty complete picture of the 
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existence and stability of localized states in canonical ID model equations, for example the Swift-Hohenberg 
model [24] with quadratic-cubic or cubic-quintic nonlinear it ies, but in contrast, comparatively little is known 
about the detailed structure of localized states in spatially discrete problems. 

Motivation for the detailed investigation of discrete models comes from a number of sources, but in 
particular discrete models arise as natural descriptions of arrays of waveguides in nonlinear optics [28 . For 
example, the discrete bistable nonlinear Schrodinger equation 

i^ n + CAljj n + S^ n \^ n \ 2 - ^ n \^ n \ 4 = 0, (1) 

for an array of complex fields i/j n (t) naturally describes properties of discrete optical 'dissipative solitons' 
in a /c-dimensional square lattice whose sites are indexed by the (multiple) suffix n G The operator A 
is an appropriate sum of nearest neighbour (and possibly next-nearest neighbour) differences depending on 
the lattice dimension coefficient k. Substituting the ansatz ifj n = u n exp(— i/it) where u n is a real stationary 
field, we see that u n must satisfy the real difference equation 

fIU n + CAu n + Sul - = 0. (2) 

Hence the possible steady states correspond to those of a discrete Allen-Cahn equation [9 , albeit with a 
cubic-quintic nonlinearity rather than the well- studied purely cubic case. 

In this Letter we investigate the existence of localized states in the canonical nonlinear lattice model ([2| 
in one and two dimensions. Such models are in many senses simpler than the corresponding spatially 
continuous equations and we discuss the similarities and differences that arise. To aid this discussion we 
add a time derivative term to ([2| which makes stability correspond to that found in the spatially-continuous 
Swift-Hohenberg model. In contrast to previous work on this canonical problem by other authors, our central 
focus in this paper is on the appearance of isolas. In particular we show that, firstly, in ID and on a finite 
lattice, isolas are an immediate consequence of changing the boundary conditions from periodic to Neumann 
(or, in fact, Dirichlet). Secondly we show that in 2D isolas exist when periodic boundary conditions are used, 
and reconnect and then disconnect again as the cell to cell coupling parameter is increased. We also comment 
on (i) the overall bifurcation structure in ID, (ii) the similarities between the ID bifurcation structure and 
that of the continuum ID Swift-Hohenberg model and (iii) the connections between ID lattice structures 
and 2D lattice structures, using the former to explain features of the bifurcation diagram we obtain for the 
latter. 

The structure of the Letter is as follows. In section [2] we introduce the ID lattice case and explore the 
bifurcation structure of localized states. We discuss the effects of different boundary conditions when the 
problem is restricted to a finite lattice. In section [3] we consider a 2D lattice, highlight the new features that 
arise, and use the ID results to explain aspects of the snaking curve. We summarise the main results of the 
paper in section [4] 

2. One-dimensional lattices 

Motivated by after applying the rescaling u n = y^u ni rescaling fi and C appropriately, dropping the 
and adding a time derivative term appropriate for strongly dissipative rather than conservative dynamics, 
we consider the following equation posed on the integer lattice Z: 

U n = C(u n+1 + U n -1 ~ 2u n ) + fjbUn + 2u\ - U*, (3) 

where u n is a real scalar variable defined at each lattice site, C is the strength of linear coupling between 
adjacent sites, and — 1 < \i < is the primary (real) bifurcation parameter. The addition of a time derivative 
term is, of course, not necessary for the investigation of the existence of equilibrium states. However, it 
is a useful addition to indicate stability and to provide a contextual similarity to the studies of the Swift- 
Hohenberg equation in spatially continuous systems referred to above. If C < then we apply the parameter 
symmetry given by applying the 'staggering' transformation 

u n ^(-l) n u n , C^-C, /i^/i-4C, (4) 
2 



(a) 



(b) 



Figure 1: (a) Bifurcation diagram in one dimension for an array of TV = 8 cells with PBC and coupling strength C = 0.2. Thick 
lines represent stable solutions, thin lines represent unstable solutions. The blue lines are asymmetric 'ladder' solutions, (b) 
The bifurcation diagram with NBC. Red curves are isolas, and blue curves terminate at branch points on the black snaking 
curve. 

which leaves the equation unchanged, but switches the sign of C. Thus in this paper we will only consider 
positive C, noting that the existence (and stability) of states {u n } with C > is equivalent to the existence 
(and stability) of states {v n } = {(— l) n u n } when C < and ji is shifted as in Q. 

Importantly, and in common with the Swift-Hohenberg equation, equation ([3| is a gradient flow, that 
is, we can write u n = —dF/du n where the potential F is given by 

F(u) = \C{u n+1 - u n ) 2 - \iiul - \ui + \ul (5) 

n 

Thus F < 0, so that every solution of ([3| flows down gradients of the potential toward an equilibrium 
solution, and every stable equilibrium state is a local minimum of the potential. No periodic or complex 
dynamics are therefore possible. 

There are five homogeneous equilibria: the trivial state u n = 0, and four non-trivial equilibria u n = ±iz±, 
where 

4 = i ± v / rT7^. (6) 

The state u n = is linearly sta ble for /i < and unstable for /i > 0. Taking the positive square roots, 
the lower uniform state U- = y/l — yl + ji exists in — 1 < /i < and is always unstable, and the upper 
uniform state u+ = y/l + yT + // exists in \i > — 1 and is always stable; there is a saddle- node bifurcation 
at // = —1. 

The per-cell potential for a homogeneous equilibrium u n = u* is given by F(u*) = —\iiu^ — + \v%. 
Since the system acts to minimize the total potential, the relative values of potential at the homogeneous 
equilibria are important. The trivial state has zero potential, and the potential of the upper state u+ depends 
only on \i. By looking for a double root of F(u) = we find that the upper state has zero potential when 
M = Mmx = — 3/4; this is the so-called Maxwell point at which the upper uniform state and the trivial state 
are energetically equal. When \i > /i mx the upper state is energetically more favourable; when \i < /i mx the 
trivial state is energetically more favourable. 

2.1. Boundary conditions 

In many physical applications it is appropriate to consider Q posed on a finite grid n G {1,...,7V} 
with boundary conditions specified with the aid of 'ghost cells' and un+i- There are three obvious 
choices: periodic boundary conditions (PBC) for which uq = un and un+i = Neumann boundary 
conditions (NBC) for which = v>i and un+i = ^at, and Dirichlet boundary conditions (DBC) for which 

^0 = ^AT+l = 0. 
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With PBC or NBC the homogeneous state is an equilibrium configuration. With DBC this is not the 
case and this obvious difference leads us to focus in this discussion only on PBC and NBC: we find there 
are substantial differences between these two cases that require some discussion. 

2.2. Symmetries 

Equation ([3| posed on the infinite integer lattice has a symmetry group generated by the independent 
operations of translation r, reflection p and an up-down symmetry n (since the nonlinearity is an odd 
function of u n ): 

(7) 

Thus the full symmetry group of the infinite lattice problem is Z k Z2 x Z2. Posing the equation on a finite 
lattice implies changes in the symmetry group. Consider the following actions of the elements r, p and n: 

T '.U n ^ U n+1 ( mod at) , p:U n ^ ^AT+l-n, k : u n i-> -u n . (8) 

With PBC the translation subgroup Z, generated by r, is simply replaced with the cyclic group Zat so that 
the full symmetry group is Z^r X Z2 x Z2 = x Z2. With NBC or DBC we lose translation symmetry 
entirely, so the symmetry group becomes Z2 x Z2 = (p) x (k). 

We will be particularly interested in localized solutions to Q which lie in Fix(p) or Fix(pr), i.e. which 
are symmetric under reflection either about a lattice point, or about a point halfway between two lattice 
points. We refer to such solutions site-centred and bond- centred. When TV is even, site-centred solutions lie 
in Fix(pr) and bond-centred solutions lie in Fix(p). When TV is odd, site-centred solutions lie in Fix(p) and 
bond-centred solutions lie in Fix(pr). 

2.3. Linear stability on a finite lattice 

(k) (k) 

The eigenfunctions Vn and eigenvalues of the discrete Laplacian, obtained by solving Avn = 
^i+i — 2i4^ + = ^ v ri^ on an TV-cell lattice with PBC are easily computed to be = e l6>fcU and 
Afc = — 4sin 2 (6 ) / c /2) where 0^ = 2irk/N for k = 0, . . . , TV — 1. Thus the zero state undergoes an instability 
to a mode with 'wavenumber' parameter Ok when p = 4C sin 2 (Ok/ '2). Since C > the first instability 
encountered when increasing p is a bifurcation to the lower uniform state U- at p = 0. This branch 
bifurcates subcritically, turning around in a saddle- node bifurcation and restabilising at p = — 1 to form the 
upper uniform state u+ as shown in figure [l] 

The upper uniform state u+ undergoes no instabilities as p is increased. In contrast, the lower state u_ 
has many secondary bifurcation points, occuring at 

P = Pfc =' ~\ ~ Csin 2 (^/2) db ^-Csin 2 (^/2), (9) 

for values of k such that the expression under the square root is positive. The bifurcations are supercritical 
'double pitchforks', each creating two new bifurcating branches of 'odd' and 'even' modulated states which 
develop into single- and multi-pulse localized states, each of which terminates on U- at both ends. In 
figure [l] only the single pulse branches are shown for clarity. This behaviour of secondary bifurcation 
creating branches of odd and even- symmetric modulated and localized states is extremely similar to that 
given for the continuum Swift-Hohenberg equation in [3j [13] without the additional complexity that arises 
from continuous changes of wavenumber along the branch which is prohibited in the present case by the 
spatial discreteness. 

Noting that the instability creating single-pulse localized states is the k = 1 branch, inspection of ([9| 
shows that the k = 1 case is the last instability to remain as C increases at fixed TV. Hence the maximum 
coupling strength above which there are no modulational instabilities of U- is C max = l/[4 sin 2 (7r/TV)]. 
Cmax therefore also provides an upper bound on the lower limit of the range of C over which localized states 
asymptotic to zero may exist. 
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Figure 2: Equilibrium solution profiles of |3j solved on a ID lattice with N = 8 using PBC and taking C = 0.2. Profiles 
correspond to the right-hand saddle-node bifurcation points of Figure Qa). The upper three solutions are site-centred and the 
lower three solutions are bond-centred. 

2.4- Homoclinic snaking 

We use the numerical continuation program auto to explore the bifurcation structure of localized so- 
lutions to Q using both PBC and NBC. Stability results come directly from auto, since ([3| is simply a 
collection of ODEs. 

With PBC the bifurcation diagram figure [I] closely resembles the homoclinic snaking picture for the 
Swift-Hohenberg equation in a finite domain, with branches of localized solutions emerging from the uniform 
branch at \i ~ —0.07 and turning around in a succession of fold bifurcations to form two intertwined snakes 
of localized structures, one of site-centred and one of bond-centred solutions, before reconnecting to the 
uniform branch at \i « —0.98. For N even, as is the case in figure [l] these solutions lie in Fix(pr) and Fix(p) 
respectively; this is shown clearly in figure [2] 

Again, just as in the continuum Swift-Hohenberg equation, the two snaking branches are connected 
by 'ladders' of asymmetric structures: the ladders comprise solutions with no symmetry at all that are 
always unstable. These branches bifurcate subcritically from the main snake, at locations near to the fold 
bifurcations making up the main body of the snake (not shown in figure [T]). Solution profiles at the lowest 
(in terms of L 2 norm) six right-hand saddle nodes are shown in figure [2J 

2.5. Neumann boundary conditions 

With PBC the bifurcation structure of the discrete model ([3| exactly parallels that of the continuous 
bistable Swift-Hohenberg equation. When we switch to NBC one of the snaking branches persists but the 
other fragments into a series of isolas (closed curves of equilibria that no longer connect up into a snake) 
and 'bubbles' (curves of equilibria that begin and end in bifurcations from the persistent snaking curve). In 
figure [I] the persistent snaking branch is indicated in black, isolas are indicated by the red curves and the 
bubbles are shown in blue. For N even, the persistent snaking curve consists of bond-centred states with 
symmetry Fix(p). The two halves of the pitchfork bifurcation where the bond-centred states bifurcate from 
the uniform state are no longer related by symmetry: in figure [ljb) we show the branch corresponding to 
states that are localized in the centre of the lattice rather than at the edges of the lattice. 

The snaking curve of site-centred states in figure [TJa) breaks up in figure [TJb) into isolas and bubbles 
which do not directly interact. In figure [ljb) the lowest isola approaches the bifurcation point where the 
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Figure 3: Regions of existence of localized states in the (/i, C) plane: 1-cell localized state (red), 2-cell (green) and 3-cell 
(magenta). The front solution is indicated by the blue curves: on an infinite lattice the edges of the snaking region asymptote 
to these curves as the width of the localized state increases. The vertical dashed line is the Maxwell point fi = —3/4 for 
reference. Results obtained for TV = 100. 



bond- centred branch connects to the uniform branch but does not bifurcate from either. As the norm 
increases, isolas alternate with bubbles which do bifurcate from the bond-centred snaking branch. The 
bubbles bifurcate from points close to where the 'ladders' bifurcate in figure [TJa) and lack symmetry: in 
fact they comprise states which evolve smoothly from the exactly bond-centred state on the snaking curve 
and become very close to the subspace Fix(pr) before moving back towards Fix(p) and reconnecting with 
the snaking branch further up. Moreover there are, properly speaking, four copies of each isola and bubble 
branch of site-centred asymmetric localized states, related by the symmetries n and p. 

This behaviour mimics what has been observed for the Swift-Hohenberg equation with non-periodic 
boundary conditions [3l[T3j[T6]. It is important to note however, that in contrast with these investigations, 
the domain size cannot be used as a continuously varying parameter in the present context. In fact, we 
anticipate that increasing TV by unity reverses the roles of the bond-centred and site-centred solutions when 
the boundary conditions are left unchanged. For example, with NBC as we consider here, when N is odd 
we expect the site-centred snaking curve to persist and the bond-centred one to fragment into isolas and 
bubbles in a manner similar to that shown in figure [I] 

3. Two-dimensional lattices 

In two dimensions our interpretation of pi must be modified since there are additional possibilities for 
the coupling terms. On the square lattice ZM;here are two natural nearest-neighbour difference operators: 

A"*" U nrn = Un-\-l,m H~ U n — i ?m -|- ^n,m+l H~ ^n,m— 1 4?i nm , (10) 

A x u nm = it n+ i 5m+ i + u n -i ifn -i + it n _i )m+ i + ii n+ i >m _i — 4:U nm . (11) 
This leads to the following 2D generalization of (|3| 

Unm ~ C + A + 1i nm + C* A X U nrn + flU nrn + 2u' nrn — V? nra ^ (12) 
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where there are now two coupling parameters C + and C x , the coefficients of nearest-neighbor (NN) and 



next-nearest-neighbor (NNN) coupling. As in the ID case, ( 12 ) is variational. For this equation our primary 
interest is in establishing that even with PBC there are isolas containing stable localized states. Therefore, 
we confine our presentation to localized solutions that are far from the lower and upper ends of the snake 
where reconnection to the uniform state occurs in a finite domain. To further simplify the discussion we 
set C x =0 for the majority of this Letter. Clearly it is of interest to investigate how the snaking structure 
evolves with C x nonzero, and we will return to this and other issues in subsequent work. 




Figure 4: Three typical localized solutions to ( |12| with C + = 0.1, C x = and /i = —0.6. The three solutions are respectively 
site-centred, bond-centred and hybrid. 

The bifurcation structure of localized solutions in 2D is made complicated by the existence of at least 
three snaking curves: a localized state can be site-centred, bond-centred or a hybrid of the two, being 
site-centred in one lattice direction and bond-centred in the other (for a more complete discussion see [6]). 
Examples of each of these types of solutions are shown in figure [4] For simplicity we will now focus on 
bond-centred states and follow the evolution of snaking curves of such states. FigureJE] shows the evolution 
of bond-centred localized states along the continuous snaking curve shown in figure |6[b). 

3.1. Bifurcation structure of localized states 

It is not at all clear theoretically how the well-established theory for homoclinic snaking in ID extends 
to 2D since there does not appear to be an analogous paradigm to the spatial dynamics approach. However, 
numerical computations show that the bifurcation diagram for localized states contains qualitative similari- 
ties (see Figure [6]). We remark that to generate this figure our computations were carried out using auto on 
a 20 x 20 grid. To accelerate the computation we simulated only a quarter of the domain, using appropriate 
'reflecting' boundary conditions to ensure that the new states generated were bond centred. This imposed 
symmetry implies that bifurcations to asymmetric 'ladder' states could not be detected. Other computations 
indicate that such bifurcations exist, and we will discuss the existence of asymmetric states elsewhere. 

The states shown in figure [5] lie on the snaking curve shown in figure [6^b) in which the bifurcation 
parameter \i is plotted on the vertical axis and the solution measure M 2 on the horizontal axis. M is a 
scaled version of the L2 norm defined as: 

/ V u 2 

Use of the solution measure M confers the advantage that at low coupling strengths it measures the (near- 
integer) number of cells that are 'active' in the localized state. Figure |6|a) indicates that 'high peaks' (e.g. 
the saddle-nodes labelled (a), (b), (c), (e) and (f) at which the value of \i is markedly less negative than 
at most saddle-nodes) are closely aligned with even integer values of M. These labels correspond to the 
labelled localized states in figure [5] that are in the form of complete 2k x 2k square arrays of active cells. 

The evolution of bond-centred localized states as we move up the snaking curve is shown in figure [5] 
Starting from[5ja) in which four cells are close to u+ and the rest are near zero, new u+ cells are added to 
the solution profile starting in the centre of the faces, and then progressing outwards towards the corners 
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Figure 5: The first 30 bond-centred localized states encountered at saddle-node bifurcation points corresponding to local 
maxima of /i, moving along the snaking curve with M increasing, for fixed C+ = 0.2. Labels (a) - (f) correspond to the labelled 
saddle-node bifurcations on the snake indicated in figure pTb). 



until the profile is square, at which point the sequence repeats itself. Numerical continuation suggests that 
it is possible to proceed arbitrarily far up the snaking curve in this manner as long as the coupling coefficient 
is sufficiently small. There are similar snaking curves of site-centred and hybrid localized states which for 
brevity, and because they are not central to our discussion of isolas, we will not discuss here. 

3.2. Snaking limits 

The horizontal lines in each part of figure [6] indicate the locations of saddle-node bifurcations far up the 
ID snake for two values of the coupling coefficient: C = C + + 2C X (red) and C = 2C + + C x (blue). These 



values are determined from the part of the blue curve in figure 3 
At small coupling strengths it is clear numerically in figure 6 
up the snaking curve, the width of the snaking region is given overall by C = C + since the existence of 
square 2k x 2k arrays of localized states is determined by the existence of stationary fronts aligned with 



which lies in \i > —0.75. 

a) that, as M 2 increases and we move far 



a lattice direction. This is easy to see from the lattice equations 12) by removing the dependence on the 
second coordinate m, i.e. setting u nm = u n we recover (|| with C = C+ + 2C X . All other saddle-node 
bifurcations involve localized states which are not complete squares and therefore involve stationary diagonal 



fronts between u+ and zero. Substituting the ansatz u nm = ug where I = n + m into ( 12 ) we obtain 



ill = 2C + (^ + i + U£-i - 2ui) + C x (ui+2 + ^-2 - 2^) + fj,U£ + 2uf - uf . (13) 

Hence, in the case C x — this reduces exactly to (|3| with C = 2C + . When C x ^ and the coupling 
coefficients are small (and hence the front is sharp compared to the lattice spacing) we anticipate that 



it£+2 ~ ~ v>+ on one side of the front and U£-2 ~ ~ on the other and so (13) can be well 

approximated by Q setting C = 2C + + C x . 

Figure [(^d) shows the bond-centred snake with C + = 0.1 and C x = 0.05. Qualitatively there is little 
difference between this and parts (a) and (b). We observe two distinct scales in the bifurcation diagram, 
with the widths of turns on the snake corresponding roughly to those in the ID model with C = C + + 2C X 
or C = 2C + + C x . We anticipate that for C x small compared with C + there will be no significant 
differences from the C x = case, but numerical continuation suggests that larger C x can alter the snaking 
curve significantly, especially when C x > C + . We should also note that with two independent coupling 
parameters there is no equivalent to the staggering transformation presented for the ID model. 

3.3. Switchbacks 

As the coupling strength C + increases the snaking curves qualitatively depart from the ID case due to 
the existence of regions in the bifurcation diagram where the snaking curve turns back on itself and the L2 
norm (or, equivalently, M) of the localized state descreases for a number of twists before turning back again 
and continuing to snake upwards. We refer to such an episode as a 'switchback'. Examples of switchbacks 
are shown in figure [7] In part (a) of figure [6] switchbacks are observed to begin only at M 2 w 64, 100 and 
144. In (b) an additional switchback appears at M 2 « 36. In (c) almost all the twists in the switchbacks 
have disappeared, although the bifurcation curves near M 2 = 64, 100 and 144 show clear similarities. 

We find that the appearance and disappearance of switchbacks as the coupling strength C + is varied 
can be explained by collisions between solution curves in the bifurcation diagram. When C + is small it is 



important to note that there are many localized solutions to ( 12 ) possessing square symmetry which do not 
appear on the snaking curve, e.g. as shown in figure [8^ a). These additional localized states exist on isolas 
in the bifurcation diagram. As C + is increased these isolas collide with the snaking branch and form a 
switchback. Assuming, as is reasonable at low coupling strength, that the isola contained a stable solution, 
the snaking curve will now contain (at least) one extra stable solution, by which we mean a localized state 
now present, but not present on the snake at lower C + . As C + is increased still further the snake narrows. 
Intriguingly, and unexpectedly, the switchback may detach itself from the snake as C + is increased further, 
re-forming a disconnected isola. Figure [8] illustrates these events with a sequence of bifurcation diagrams 
in the (/i, L2) plane. The processes of attachment and detachment are facilitated by cusp singularities in 
which pairs of saddle-nodes on the snake collide and disappear as the coupling strength increases. Figure [9] 
shows the cusp bifurcations in the (/i,C + ) plane corresponding to the attachment and detachment of the 
isola near M 2 = 36 shown in figure [?^a) and figure [SJ 

In summary, these isolas provide a mechanism by which additional turns can be added to the snake, 
and they explain the unintuitive switchback episodes of the snaking curve as the localized states increase 
in size. This process of isola attachment and detachment typically occurs over a remarkably small range of 
C + , indicated by the vertical scales in figure [9J At larger L2 norms the complete bifurcation diagram in the 
(/i, C + ) plane is complicated since many isolas attach and detach from the snake at nearby values of C + . 



4. Summary 

In this Letter we have discussed the bifurcation structure of a model lattice equation, the discrete Allen- 
Cahn equation with cubic and quintic nonlinear it ies, in one and two dimensions. This equation is of interest 
in its own right as a model for spatially discrete pattern formation as well as being closely related to coupled 
nonlinear Schrodinger equations which have received substantial recent attention in nonlinear optics. A 
subcritical instability of the trivial state, and the resulting bistability, leads to the creation of localized 
states which are stable over an open interval of the parameter /i, stabilized by being pinned to the lattice. 
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As remarked on by previous authors, these localized states arrange themselves into homoclinic snaking 
curves. In a finite domain we demonstrated the existence of a maximum coupling strength above which 
there is no modulational secondary instability which would allow creation of localized states: in systems 
with a sufficiently large coupling strength localized states do not appear. 

In ID we showed that breaking discrete translational symmetry, for example by using Neumann boundary 
conditions, allows one of the snakes to persist, but cuases the other to fragment into a sequence of isolas and 
bubbles ('S'-shaped curves bifurcating from the persistent snake). This contrasts with the situation in 2D 
where isolas exist even when discrete translational symmetry is not broken. These isolas attach and detach 
themselves from the primary snaking curve as the linear coupling strength is varied, creating kinks in the 
snaking curve, where the solution norm decreases along a part of the curve before beginning to increase 
again, that we refer to as 'switchbacks'. 

Numerical investigations suggest quantitative relations between the widths of the snaking curves in ID 
and 2D. For wide localized states in 2D the values of fi that limit the region over which localized states 
exist correspond to the limiting values for fronts aligned with lattice directions or on the diagonal. This 
is the natural interpretation and agrees with the observations and computations by Lloyd et al. [20] for 
localized patches of hexagons in the 2D Swift-Hohenberg equation. This correspondance between ID and 2D 
is asymptotically exact in the limit of small coupling coefficients: higher-order corrections arise at non-zero 
coupling and from the corners of the wide localized state where the fronts meet. 

Our results are closely related to the work of Carretero- Gonzales and Chong et al. [8] on the discrete 
nonlinear Schrodinger equation with cubic and quint ic nonlinear it ies. They discussed the time- independent 
Allen-Cahn equation from a spatial dynamics viewpoint, explaining the one-dimensional snaking curve as 
arising from intersections of the stable and unstable manifolds of a hyperbolic fixed point in a two-dimensional 
map. They also explored a variational approach, minimizing an effective Lagrangian to give approximate 
localized states in one and two dimensions. 

Our work considers additional mechanisms that complicate the usual snaking picture by causing isolas 
to appear: these include the effect of boundary conditions in ID and exploring larger localized states, higher 
up the snaking curve in 2D where new effects, including the switchbacks, arise. 
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Figure 6: Snaking of bond-centred solutions to jl2} with (a) C+ = 0.1, (b) C+ = 0.2, (c) C+ = 0.3 (in all cases C x = 0) and 
(d) C+ = 0.1, C x = 0.05. Note the switchbacks in (b) occur near M 2 = 36, 64, 100 and 144, see also Figure^ Solid black lines 
indicate the Maxwell point /i = -0.75. Red horizontal dashed lines at (a) \i w -0.41925, (b) \i w -0.56076, (c) fi w -0.63832 
and (d) /i w —0.56076 indicate the asymptotic location of saddle-node bifurcations on the ID snake with coupling strength 
C = C+ + 2C X , far up the ID snake. Blue horizontal dashed lines at (a) fi w -0.56076, (b) /i « -0.68309, (c) /i w -0.72479 
and (d) /i w —0.60495 indicate the asymptotic location of saddle-node bifurcations on the ID snake with coupling strength 
C = 2C+ +C X . Thin and thick lines indicate unstable and stable solutions, respectively. Labels (a) - (f) in part (b) correspond 
to the solutions labelled in figure [5] 
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(c) C+ = 0.27 (d) C+ = 0.28 

Figure 8: (a) The snaking curve (black) and an unconnected isola (red) at C + = 0.19, (b) topology change in the snaking 
curve as we increase the coupling strength, creating a switchback, (c) narrowing of the snaking curve as coupling strength is 
increased further, (d) eventual detachment of the switchback to form another isola. 
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-0.716 



Figure 9: The process of isola attachment and detachment in the (/i, C + ) plane by following limit points: on the left an isola 
attaches itself to the snake as C + is increased, briefly creating a situation with three limit points, two of which disappear in a 
cusp singularity [corresponding to Fig . |8[a)— (b)], and on the right the same process in reverse, as an isola detaches at a higher 
value of C + [corresponding to Fig. |8rcj"^(d)]. 
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